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Abstract 

The neuromagnetic activity (magneto-encephalogram, MEG) from healthy human brain and from 
an epileptic patient against chromatic flickering stimuli has been earlier analyzed on the basis of a 



memory functions formalism (MFF). Information measures of memory as well as relaxation param 



eters revealed high individuality and unique features in the neuromagnetic brain responses of each 
subject. The current paper demonstrates new capabilities of MFF by studying cross-correlations 
between MEG signals obtained from multiple and distant brain regions. It is shown that the MEG 
signals of healthy subjects are characterized by well-defined effects of frequency synchronization 
and at the same time by the domination of low- frequency processes. On the contrary, the MEG 
of a patient is characterized by a sharp abnormality of frequency synchronization, and also by 
prevalence of high-frequency quasi-periodic processes. Modification of synchronization effects and 
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dynamics of cross-correlations offer a promising method of detecting pathological abnormalities in 
brain responses. 
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I. INTRODUCTION. SYNCHRONIZATION AND COLLECTIVE EFFECTS IN 
TIME SERIES ANALYSIS OF COMPLEX SYSTEMS 



One of the main factors determining the evolution of complex systems is the presence 
of collective effects arising from an interacting or redistributing of the certain connections 
between parts of composite system. In many cases it is impossible to make an adequate 
analysis of the functions of the systems by ignoring the underlying collaborative mechanisms. 

There are various approaches used in studying the collective phenomena in complex sys- 
tems. Somehow or other, all of them are based on the analyzing unique features of the 
connected systems: certain quantitative and qualitative ratios between the system elements, 
a dynamic coordination of components under the external influences, specific synchroniza- 
tion phenomena. Some recent results have been derived by studying the effects of frequency 
and phase synchronization jl-5]. These methods are based on revealing the characteristic 
frequencies and analyzing the differences in the phases of dynamic variables derived by means 
of the Fourier transform, Hilbert transform 1, 2, and wavelet-transform [5]. Within the 
framework of another methodology the stochastic synchronization is studied by comparing 
topological structures of attractors, describing the dynamics of two nonlinear coupled oscilla- 
tors [6]. The "generalized synchronization" relationship {t| also uses the topological method 
and is the successful original development of the stochastic synchronization approach. 

Another approach to study the collective effects in complex systems is the analysis of 
cross-correlations, i.e. the probabilistic relation between the sequences of random variables. 
The cross-correlation method is used to describe the collective phenomena in various systems 
(physical, economical, biological and physiological). The perspective approach in this field 
of research is detrended cross-correlation analysis, which has been introduced in [8j] to study 
the power-law cross-correlations between nonstationary time series of various nature. This 



method is widely used in financial systems |9l-|ll|. 



There are several alternative methods to analyze cross-correlations, particularly, random 

n n n n 

matrix theory [12|, approximate entropy [13] and sample entropy [14j. In [12], the authors 
used random matrix theory, a method originally developed to study the spectra of complex 



nuclei, to analyze the mutual dynamics in price changes of the stocks. In 



131 ]. the authors 



used the method of approximate entropy, a model independent measure of sequential irreg 



ularity which is based on Kolmogorov entropy, as an indicator of system stability. In 



the authors proposed the sample entropy, a modified and unbiased version of approximate 
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entropy, as a measure of degree of asynchrony in physiological signals. 

Although these earlier methods find useful applications in real-life time series, there 
exists another set of methods in cross-correlation analysis which uses the cross-correlation 
functions itself (3), [l6|- In paper the authors compared the linear synchronization 
measures including the cross-correlation functions to nonlinear ones and revealed that for 
the considered experimental data all measures ranked the synchronization levels of the three 



examples in the same way In work 



15], the authors develop a theory of neuronal cross- 



correlation functions for analyzing the neuron interactions in large neural networks including 



several highly connected sub-populations of neurons. In work [16j, the authors developed a 
method based on the model of fractional Brownian motion and Hurst exponent to describe 
the coupling of the non-stationary signals with long-range correlations. 

In this study, we demonstrate the fundamentally new opportunities of studying the col- 
lective effects in complex systems. Our method is based on generalization of the memory 



functions formalism 



17 



18j in a case of cross-correlations between the spaced elements of 
the studied system. The important advantage of this method is the description of the cross- 
correlations in different relaxation scales in the time series of a complex systems. Here, we 
consider neuromagentic responses (magnetoencephalogram, MEG) of brain as a suitable ex- 
ample of a time series of a complex system, i.e. the human brain. Earlier analysis of MEG 
signals jl9] were performed on the basis of the memory functions formalism (MFF), and 
revealed an important autocorrelation difference between MEG signals of healthy subjects 
and of a patient with photosensitive epilepsy (PSE). Particularly, this difference appears in 
qualitative alterations of the power spectra of memory functions. Besides, it has been shown 
that the statistical memory effects play a key role in identification of PSE. 

Here we investigate cross-correlations in the MEG responses simultaneously obtained 
from multiple brain regions. We will show that mechanisms of formation of the PSE are 
connected, first of all, with abnormality of interrelations between the spaced areas of a 
cerebral cortex, which result in suppression of its regulator functions at formation of the 
response to external influences. The paper is structured as follows. Section 2 presents the 
basic relations of the memory functions formalism for cross-correlations. Section 3 details 
the experimental details. Section 4 contains our results including the calculation of cross- 
correlation functions, memory functions and their power spectra. Section 5 offers the general 
conclusions about connections between the PSE pathological alterations and suppression of 
coordination effects (frequency-phase synchronization). 



II. BASIC RELATIONS OF THE MEMORY FUNCTIONS FORMALISM IN 
CASE OF CROSS-CORRELATIONS 



Following 17H20|. we consider stochastic dynamics of the magnetic induction gradient, 
registered in two different brain regions as the sequences {xj}, {yj} of random values X, Y: 



(1) 



X ={x(T),x(T + t),x(T + 2t),..., x(T + (N - l)r)}, 
Y ={y(T),y(T + r),y(T + 2r),..., y(T + (N - l)r)}, 

where T is the initial time point, (N — l)r is the time period of signal registration, r is the 
time interval of signal discretisation. 

Mean values, fluctuations and dispersions for a set of random values (CQ) can be written 
as follows 

^ N-l „ N-l 

( X ) = N S X ( T + JT ^ X i = X ^ T + JT )' Sx i = X i ~ a * 

3=0 

N-l N-l 

< y >= «E»( r + J T )' Vi = v(T + jr), Sy, = yi -(Y), oj = - £ 6$. 



j=0 



3=0 j=0 

To describe the probabilistic relation between the sequences of random variables X and 
Y we use the normalized time-dependent cross-correlation function (CCF): 

N-m-l 



C(t) 



1 



. — y\ Sx(T + jr)8y(T+(j + m)r), 
(N - m)a x a y 



(2) 



t = rnr, 1 < m < N — 1. 
Function fl2l) satisfies the conditions of normalization and relaxation of correlations: 



Iimc(t) = l, limc(t) = 0. 

t— ¥0 t— >00 

It should be noted that the second property not always satisfied for the real time series of 
complex systems. 

Taking into account the fact that the process is discrete, we must rearrange all standard 
operation of differentiation and integration 2jJ, |22| : 



dx Ax(t) x{t + r) — x(t) 
~dt ~^ At ~ t 



/l) n—1 n—1 

x(t)dt = x{T a + jr)At = r x{T a + jr) = nr(X), b — a = c,c = n.r 

j=0 j=0 



(3) 



Following 17H19|. to describe the dynamics of studied variables let's introduce the evo 
lution operator U(t + r, t): 

x(t + r) = U(t + r,t)x(t), 
yit + r) = U(t + r,t)y(t). 
Using Eqs. ([3]) and (Tjl]), we can establish a formal equation of motion for any Xj and yf 
Axj (t) Xj + i (t + r) — Xj (t) 



(4) 



At T 

A yj(t) = yj+i(t + r) -yj(t) 

At T 



T- l {U(t + T,t) - l}Xj(t), 



(5) 



{U(t + r,t)-l} yj (t) 



and then introduce a Liouville's quasioperator L as follows: 

dx(t) Ax{t) 



-'r ; - iL(t,T)x(t), 

dt ^^T = ^(t,r)y(t), 
L(t,r) = (tr)- 1 [U(t + T,t)-l]. 



dt At 
dy(t) ^ Ay(t) 



(6) 



Analogously by [17H20|. let's present the sets of fluctuations 5xj = 5x(T + jr), 5yj = 
5y(T+jr), where j = 0, 1, . . . , N— 1, as the A;-component vectors of system state in Euclidean 
space of the vectors of state: 



A° = A°(0) = {Sx , 5xi, ... , fafc-i}, 
B° k = Bl(0) = {8y ,8y 1 ,...,8y^ 1 }. 



(7) 



^-component vectors of system state in time moment t: A™ +k (t), B™ +k (t) are arose by the 
time shift on step t = mr: 



A m+fc — A m+k {t) — {5x m , 5x m +i, . . . , 5x m +k-i} i 



(8) 



B m+ fc = B m+fc (*) ={Sy m , 5y m+1 , . . . , <Wfc-i}- 

These vectors can be derived by multifold actions of the evolution operator U(t + r, t) to 
vectors of the initial system state: 



A™ +k (t) = U(T + mT,T)Al(0), 



B™ +k (t) = U(T + mT,T)Bl(0). 



(9a) 
(9b) 



Further we will introduce the scalar product operation: 



8=1 



Using the Eqs. ([2]) - (|9b|) . we can rewrite equation for CCF as a scalar product of vectors 
A°(0) and B- +fe (t): 

m _ (A°(0) U(T + mr,T) B°(0)) _ (A°(0) B^ fc (t)) 

cl j (A2(o) B° fc (o)> (Ag(o) Bj(o)> • 1 Uj 

Let's write the finite-difference Liouville's equation for the state vectors (see (jSJ): 

^B™ +fc (t) = zL(t,r)B™ +fe (t). (11) 



Let's introduce the operators of projection II and P for the Euclidean space of state 
vectors: 

n = |B°(0))(A°(0)| 2 = = _ 2 = = = 

(A2(o) b2(o)> 1 } 

The initial CCF ({TO]) can be derived by projection of state vector B™ +fc (t) to vector of 
the initial system state B° k (0): 

nB- +fc (t)> = B °(o)>^|^|^ = B° fe (0)>c(t). 

Projection operators II and P split the Euclidean vector space into two mutually-orthogonal 
subspaces: 

B(fc) = B'(fc) + B"(k), B'(k) = nB(fc), B"(fc) = PB(fc), B™ +k (t) e B(k). 

This permits to split Liouville's equation f lTTj) into two equations within two mutually sup- 
plementary subspaces as follows: 

AB'(t) 



At 
AB"(t) 



iL n B'(t)+iL 12 B"(t), (13a) 
iL 21 B\t) + iL 22 B"{t). (13b) 



At 

Here = UiLUj, i,j = 1,2, ITi = II, II2 = P, are the matrix elements of Liouville's 
quasioperator: 

L = Ln + h\ 2 + L 2 \ + Z/22- 

To transit between the subspaces the operators are used in the following way: 

L11— from B' to B', Li 2 — from B" to B', L 21 — from B' to B", L 22 — from B" to B". 

Solving the Eq. ( 113b|) and using the derived results in Eq. (I13a[) we come to the closed 
finite-difference discrete equation for the initial CCF (se for details [l3]): 

^ = Afc(t) -rArj^MrUrW-jr). (14) 

j=0 
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Here Xf Y is an eigen frequency of Liouville's quasioperator L, A^ Y is a relaxation parameter 
with the dimension of squared frequency, M^ Y (jr) is a normalized cross-correlation memory 
function of the first order: 

x xy = (A°(0)LB°(0)) XY (A°(0)L 2 B°(0)) 
(AJ(O)BJ(O)) ' 1 (AJ(0)B2(0)> ' 

1 U ] (A2(0)BJ(0)> 1 U 1 J 

Using Eqs. tTTOj) - ([13]) . we can introduce a set of projection operators Il n and P n and 

derive a chain of equations for cross-correlation memory functions of the n — 1 order: 

= ^ Y M^(t) - rA;T £ M* Y Ur)Ml\(t - jr). (16) 

j=0 

To derive the kinetic and relaxation parameters and also memory functions we will use 
the Gram-Schmidt orthogonalization procedure: (W^W y ) = 5 ntm (W^W Y ), where 8 n>m 
is Kronecker's symbol. Using it we may easily introduce the recurrence formula for the 
dynamic orthogonal variables Wjf , W y in which the senior values are connected with the 
junior values: 

W* = A°(0), W? = (iL - Af)W^, = (it - Af)Wf - Af - . . . , 
W y = B°(0), Wr = (it - Xf Y )W Y , W Y = (it - \* Y )W Y - Af Y W Y . 

Then the eigenvalues of Liouville's quasioperator \* Y and relaxation parameters A* Y in 
Eq. (USD will be: 



Normalized cross-correlation memory function in ([TBI will be: 

n ~ n ; (W x W y ) 

\ n n I 

A relaxation time of the initial CCF and memory functions of the n order are determined 
as follows: 

N-l N-l 

r c = At ]T <t,), ... , t m xy = AtJ2 M^it,). (17) 

3=0 3=0 

The set of dimensionless values will determine the statistical spectrum of non-Markovian 
parameter, the informational measure of memory: 

( r XY\_fXY c XY r XY\ e XY _ T c XY _ Tm *-i /-, o\ 

T M? Y T MXY 



Thus the value s^_ x is an useful criterium for comparison of relaxation times of the memory 
functions M^2[ and M* Y . This quantitative criterion characterizes quantitatively the degree 
of Markovity of the processes and the memory effects in discrete dynamics of complex 
systems. 

In this paper we use the frequency-dependent case of information memory measure: 



.XY 



(19) 



Here fif (u) is a power spectrum of zth memory function: 



N-1 



cos lnvt, 

3=0 



N-1 



N-1 



At ^2 M i Y (tj) COS 2nut 3 

j=0 



At M * Y {tj) cos 2-Kvtj 

j=0 

The preliminary analysis of experimental series of complex systems shows, that the special 
interest for analyzing the statistical memory effects represents the ultralow frequency area 
of the information characteristic e^iy — > 0): 



^P(o) 



(20) 



This area determines the long range correlations in experimental series {xj}, {i/j}- 



III. REGISTRATION OF NEUROMAGNETIC ACTIVITY IN HUMAN CERE- 
BRAL CORTEX 



The analyzed MEG signals 23|, |24| are the induced neuromagnetic responses obtained 
from a group of neurologically healthy participants and of a patient diagnosed with PSE 
while they were viewing flickering stimuli with various chromatic combintions. The PSE is 
a form of reflexive epilepsy in which the seizures are provoked by various forms of visual 
stimuli. Amongst various parameters of a visual stimulus, chromaticity is less studied in 
the context of PSE, yet flickering colourful stimuli are quit e widespread in modern technical 



age laden with multimedia gadgets. In previous works |23l. |24J| . the authors used the whole- 
head MEG system (Neuromag-122, Neuromag Ltd. Finland) by means of 61 SQUID-sensors 
(superconducting quantum interference device) and recorded neuromagnetic brain responses 
from nine healthy or control subjects (age range 22-27 yrs) and a patient with PSE (age 12 



yr). Subjects of control group had no personal or family history of photosensitive epilepsy 
All subjects were explicitly informed that flicker stimulation might lead to epileptic seizures. 
In order to limit the health risk, stimuli with only short duration (2 s) was used. All 
participants gave their written informed consent before recording. Adequate clinical protocol 
was maintained during the recording from the patient. The subjects were instructed to 
passively observe visual stimuli with minimal eye movement. Visual stimuli (red-blue, red- 
green) were generated by using two video projectors (Sharp XV-E500, Japan), each of which 
produced a continuous single color. For each trial, the flickering stimulus was presented for 
2 s; the gap between trials was 3 s. Evoked MEG responses were obtained by traditional 
averaging technique across artifact-free trials; at least 80 artifact-free trials were averaged 
for each stimulus. The sampling frequency was 500 Hz. 

To study cross- correlations from various combinations of the MEG signals, we have chosen 
a subset of sensors from selected brain regions as follows: occipital (sensor no. 51, 52, 53), 
left-temporal (No. 30, 34) and right-temporal (No. 56, 57). Our choice is determined by the 
physiological mechanism of processing and transferring the visual information in a brain. 
The visual cortex, in which the signals are transferred from the retina via thalamus, is in 
occipital area. Temporal areas play an important role in the perception of the information, 
including visual one. Out of three chromatic flickering stimuli (Red-Blue, Blue-Green, and 
Red-Green), we have chosen to analyze the MEG signals against Red-Blue flickering stimulus 
only, as this was shown to cause the largest neuromagnetic responses {23, 24]. For the 
generalized description of dynamics of the cross-correlations we consider the sixth subject 
from control group. 



IV. CROSS-CORRELATIONS IN INDUCED NEUROMAGNETIC ACTIVITY OF 
HUMAN BRAIN 

A. The PSE-induced stratification of phase clouds of the dynamic orthogonal 
variables 

Figs. 1, 2 contain the initial time series registered by the 56th (right-temporal area) 
and the 51st (occipital area) SQUIDs of the healthy subject (Fig. 1) and the patient with 
PSE (Fig. 2). MEG dynamics of healthy subject is characterized by the significant large- 
scale fluctuations while in patient the small-scale fluctuations are manifested against the 
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background of quasi-periodic oscillations. Quasi-periodic structure of the signal in Fig. 1 
is directly connected with the physiological rhythms of electromagnetic brain activity in 
healthy subject. In the first 200 ms the flickering stimulus is switched off, the control signal 
is registered. It is characterized by smaller amplitude of the fluctuations in initial time series. 
Switching-on of the flickering stimulus causes the increasing of average value of signals in 
healthy subject, and in patient's signals that causes the more significant fluctuations. 

Figs. 3, 4 demonstrate the plain phase projections of the phase clouds of dynamic or- 
thogonal variables {W^, Wf}, where i = 0...3. To represent the generalized view of phase 
portraits of possible combinations of the dynamic orthogonal variables {Wj^ 3 , 3 } we 
have chosen first four of them. On phase portraits for healthy subject's brain signals (Fig. 3) 
two phase areas similar on structure are distinctly manifested. Smaller of them (it is shown 
by the arrow) corresponds to the magnetic field dynamics when the flickering stimulus is 
switched off. Switching-on of stimulus causes the occurrence of transition between the phase 
cloud areas that is connected with changes of amplitude in initial signals (Fig. 1). 

Other picture is observed in the phase portraits of brain signals in patient with PSE 
(Fig. 4). First, phase clouds have bigger sizes than those for the healthy subject. Secondly, 
they have stratified structure concerning the central nucleus. It means, that switching-on of 
stimulus does not cause the appreciable phase changes in dynamics of the brain signals. 

B. Spectral cross-correlation properties of human neuromagnetic responses 

Power spectra of used cross-correlation functions for the studied neuromagnetic signals 
allow to extract unique information about the collective phenomena in brain functioning. 

To demonstrate the general mechanisms of coordination in cerebral cortex signals we 
present power spectra fj,* \v) , $ (u) of auto correlation functions for two signals and power 
spectrum of respective cross-correlation function Hq Y (u) in Fig. 5. These dependencies 
are derived for neuromagnetic responses of one of the healthy subjects. Their comparison 
distinctly detects the increasing of respective peaks in CCF power spectrum (see, for exam- 
ple, peaks 1 and 2). Actually, we reveal here a frequency-phase synchronization between 
the cerebral cortex signals as a response to color flickering influence. Thus analysis of the 
initial CCF power spectra allows to detect synchronization of the brain signals and reveal 
its frequencies. 

Power spectra of initial CCF and respective memory functions for the brain signals of 
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healthy subject (Fig. 6) are characterized by two peaks. The most significant of them is 
in the frequency 9.18 Hz, i.e. it is within the frequency range of normal brain rhythms. 
These rhythms reflect the complex psycho-physiological processes of brain activity. There- 
fore strong change of their characteristic frequencies is the indicator of functioning abnor- 
malities in human brain and the central nervous system. The second peak has frequency of 
17.56 Hz. Besides memory functions power spectra contain the additional peaks in higher 
frequencies. They reflect specific features of the brain bioelectric activity. Note, that their 
frequencies are multiple of 8.3 Hz, i.e. of frequencies of the normal brain physiological 
rhythms. The analysis of CCF and memory functions power spectra derived for combina- 
tions of control group signals, shows, that all principal peaks are in the field of low frequencies 
(up to 50 Hz), and the highest peaks have frequency of 8 - 12 Hz. Thus, the frequency-phase 
synchronization of brain signals in healthy people occurs in a low-frequency range. 

Just the opposite picture is observed in power spectra of cross-correlation functions for 
brain signals of the patient with PSE (Fig. 7). Peaks fill in all range of frequencies. Thus 
the process with frequency of 49.33 Hz is dominating. Spectra in all relaxation levels have 
the same kind. 

Specific character in interaction of brain signals of the patient with photosensitive epilepsy 
reflects abnormalities in reaction to flickering stimulus. Considered power spectra allow to 
take apart the development of abnormal high neuron collective activity of spaced cerebral 
cortex pathological reaction to visual influences, resulting in epileptic seizure. 

At the same time in healthy people we reveal the original protective mechanism blocking 
the development of such reaction. At PSE the suppression of this mechanism is reflected 
in domination of more high-frequency processes on intensity comparable with the normal 
physiological rhythms. 



C. Differentiation of intensity of the statistical memory at PSE 



To reveal the crucial role of statistical memory effects in cerebral cortex activity at PSE 
we shall use the cross-correlation information measure of memory ef Y (0). At first the non- 
Markovian parameter e was introduced in statistical physics of condensed matter 25|, [26] , 
however later it proved itself in detecting the physical mechanisms of abnormal functioning 
of live systems 17H19|. 

Both non-Markovian parameter e and its generalization e^- Y (0) for cross-correlation anal- 



11 



ysis allow to characterize the degree of statistical memory effects in long-range component 
of discrete dynamics. If y (0) ^> 1, then the Markovian components dominate in dynamics 
of stochastic processes. In this case the time of memory existence is too much shorter than 
relaxation time of initial CCF. The weak (short) statistical memory is manifested. Decrease 
of this information measure characterizes the expansion of memory lifetime. e^ y (0) ~ 1 
means that the processes are characterized by long-range (strong) statistical memory. In 
this case the time of memory existence is covariant with the relaxation time of initial CCF. 
When ef y (0) > 1 we can consider the studied processes to be quasi-Markovian with moder- 
ate (intermediate on time existence) statistical memory. Thus, the introduced quantitative 
criterium allows to parameterize the memory effects intensity and the velocity of loss of 
relaxation processes. 

In Figs. 8, 9 we have presented the frequency spectra of information measures of memory 
e 5i-56^^ w h ere i — i...3 ; f or neuromagnetic responses of spaced SQUIDs No. 51, No. 56. 
Statistical memory intensity in time intervals about length of initial series is determined by 
the low-frequency areas, and the high-range frequencies correspond to smaller time scales. 
To study the general manifestations of statistical memory effects in initial time series we shall 
consider only the parameter e^" 56 (0). For healthy subject the value of this parameter is 9.15. 
It means that the statistical memory effects are decreased in scales of time series length at 
the same time. For patient the parameter e^ 1_56 (0) value is 1.61. Frequency dependences of 
s^ m (v) in case of healthy subject (see Figs. 8b, c) are characterized by the significant peaks 
in low- and high-frequency ranges against the background of general oscillation structure 
with aliquot frequencies. At PSE (Figs. 9b, c) such structure of information measures 
spectra e^ -56 ^) is depressed. 

To study the general character of statistical memory manifestation in the dynamics of 
neuromagnetic responses of spaced cerebral cortex areas we derived the ratios of averaged 
values of the information measure ef Y (Q) for control group and respective value for pa- 
tient with PSE. These ratios are derived for all combinations of considered SQUIDs and 
are presented in Table 1. Analysis of the presented values allows to divide the SQUIDs 
combinations into three groups: 

1. SQUID combinations: 52-30, 52-34, 52-56, 52-57. The differences in information mea- 
sures ef Y (0) are less than of equal to 2.5 times. At PSE the coupling of signals, 
generated by these cerebral cortex areas, considerably does not change. 
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Table I: Ratio of averaged value of information measure £^^(0) for control group and respective 
value for patient with PSE. Considered SQUIDs are located in occipital (No. 51, 52, 53), left- 
temporal (No 30, 34) and right-temporal (No. 56, 57) areas of the cerebral cortex 



SQUID 


30 


34 


56 


57 


51 


4.22 


24.10 


6.03 


20.48 


52 


0.92 


1.45 


0.72 


2.42 


53 


3.69 


37.35 


4.28 


3.40 



2. SQUID combinations: 51-30, 51-56, 53-30, 53-56, 53-57. At PSE the mutual dynamics 
of neuromagnetic responses, generated by these areas is characterized by the significant 
increase of memory effects, that may be considered as the diagnostic indicator of 
pathological changes. 

3. Finally the mutual behavior of signals, registered by SQUIDS 51-34, 51-57, 53-34 is 
characterized by the most dramatic (in 20 - 37.5 times) increase of statistical memory 
intensity at PSE. Coupling of these areas changes most strongly. 

Actually, parameter ef Y (0) allows to quantitative differ the manifestation of statistical 
memory effects in mutual dynamics of neuromagnetic responses in patient with PSE. Thus 
the possibility of revealing the peculiar zones of anomalous collective neuron excitation (as 
an example, SQUID No. 51) at PSE, at which the depression of brain regulator functions, as 



a response to nickering color stimuli occurs. Note that earlier in paper 19j some conclusions 
have been made about the PSE's abnormalities of MEG signal co-ordinations at influence 
of flickering stimuli. 



V. CONCLUSION. THE ROLE OF CROSS-CORRELATIONS AND COLLEC- 
TIVE EFFECTS IN COMPLEX SYSTEMS 

Time signals generated by a complex system contain the unique information about its 
organization and character of couplings between its constituent systems. Therefore a repre- 
sentation of this information as a set of variables, i.e. parametrization of complex systems 
dynamics, is crucial to understand the complex system. To analyze the quantitative and 
qualitative parameters of the mutual dynamics of complex system, we here use the cross- 
correlation functions. 
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To study the mechanisms of collective phenomena in dynamics of cerebral cortex, we 
analyze MEG signals from healthy healthy subjects and a patient with PSE while viewing 
chromatic flickering stimuli, we have used generalization of memory functions formalism in 
case of cross-correlation analysis. 

Analysis of cross-correlations and collective effects in the mutual dynamics of MEG re- 
sponses reveals a protective mechanism present in healthy brain, which acts against adverse 
external perturbations and manifests in certain regime of frequency-phase synchronization. 
In healthy subjects this mechanism blokes the developing of epileptic reaction to the influ- 
ence of flickering light. In this case power spectra of the initial cross-correlation function 
and memory functions demonstrate a clear dominance of the low-frequency processes in 
coupling of MEG responses, respective to the normal brain rhythms. At the same time 
the high-frequency processes become dominant in the MEG responses brain signals of the 
patient with PSE. Therefore it is reputed that the substitution of low-frequency synchro- 
nization of the MEG signals by high-frequency synchronization at flickering light influence 
is an reliable indicator of PSE. 

The proposed analysis has offered not only a new insight into the character of the patho- 
logical changes in patient's brain responses but also reveals a high degree of individuality in 
neuromagnetic brain responses. 
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Figure 1: Induced neuromagnetic responses to the red-blue nickering stimulus from the occipital 
(51-st SQUID) and the right-temporal (56-th SQUID) cerebral cortex areas of healthy subject. 
The large-scale fluctuations against the quasi-periodic background are connected with manifesting 
the physiological rhythms of brain bioelectric activity 
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Figure 2: Magnetoelectric signals, generated by occipital (51-st SQUID) and right-temporal (56- 
th SQUID) cerebral cortex areas of patient with PSE under the influence of red-blue flickering 
stimulus. Signal structure is determined by the superposition of the small-scale fluctuations against 
the quasi-oscillatory background 
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Figure 3: Typical structure of plane portraits of the dynamic orthogonal variables Wf = /(Wjf ), 
where % = 0...3, for MEG signals of healthy subject. Areas of phase clouds, indicated by the arrow, 
correspond to neuromagnetic responses in absence of external flickering stimulus 
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Figure 4: Plane projections of phase portraits of the dynamic orthogonal variables {Wff , Wf }, 
where i = 0...3, for the induced cerebral cortex signals from patient with photosensitive epilepsy. 
The apparent increasing of scales for phase clouds is revealed 
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Figure 5: Power spectra of auto-correlation functions for signals X and Y in comparison with 
power spectra of cross-correlation function c(t). Frequency synchronization of these signals results 
in increasing of respective peaks on power spectrum of cross-correlation function. 
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Figure 6: Typical view of power spectra of the initial cross correlation function and memory 
functions for the mutual dynamics of neuromagnetic responses of healthy subject. Peaks, which 
characterize the periodic features of brain neuromagnetic activity, are in a low-frequency range of 
physiological rhythms 
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Figure 7: Power spectra of the initial cross correlation function and memory functions for induced 
MEG signals from patient with photosensitive epilepsy. As it is seen from the specific low-frequency 
dynamics in healthy subject, the suppression of normal physiological rhythms by high frequency 
processes is observed in this case. The frequency synchronization abnormality can be interpreted 
as a diagnostic indicator of PSE pathological manifestation. 
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Figure 8: Frequency dependences of the information measures of memory ef (u), where i = 
1...3, for signals, registered from 51-st and 56-th SQUIDs. Mutual dynamics of these signals is 
characterized by quasi-Markovity and sufficient weak intensity of statistical memory. Oscillatory 
type of spectra of the information measure £2,3^ (f ) is connected with manifesting the brain rhythms 
in MEG signals of healthy subject 
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Figure 9: Frequency spectra ef Y {y), where i = 1...3, for neuromagnetic responses from cerebral 
cortex of patient with PSE. The typical alteration of the intensity of statistical memory in dynamics 
of patient's neuromagnetic responses to external light influence in comparison with control group is 
observed. According to this fact the effects of statistical memory play a crucial role in pathological 
changes in case of PSE 
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